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Abstract 

We study collective excitations in systems described by chiral kinetic theory in external magnetic 
field. We consider high-temperature weak-coupling plasma, as well as high-density Landau Fermi 
liquid with interaction not restricted to be weak. We show that chiral magnetic wave (CMW) 
emerges in hydrodynamic regime (at frequencies smaller than collision relaxation rate) and the 
CMW velocity is determined by thermodynamic properties only. We find that in a plasma of 
opposite chiralities, at frequencies smaller than the chirality-flipping rate, the CMW excitation 
turns into a vector-like diffusion mode. In the interacting Fermi liquid, the CMW turns into the 
Landau zero sound mode in the high-frequency collisionless regime. 


1 



I. INTRODUCTION 


Parity-violating transport in the plasma of chiral fermions originating from chiral anomaly 
and driven by magnetic field attracted significant amount of interest recently. Such phe¬ 
nomena can be realized, for example, in the early stages of heavy-ion collisions at RHIC and 
LHC where a very strong magnetic field is created by large ultrarelativistic nuclei [I]. One 
such effect which may have a significant phenomenological importance is the Chiral Mag¬ 
netic Effect (CME) [IHS]. The CME causes vector (axial) current along the direction of the 
magnetic field in the presence of axial (vector) chemical potential. In heavy ion collisions 
this should lead to dipolar charge separation perpendicular to the reaction plane, whose 
experimental signature may be observed in two particle azimuthal correlations of charged 
particles IMD]. The CME has been demonstrated in various theoretical frameworks, such 
as in hydrodynamics EH, in lattice QCD [T2HT6] . in effective theories mm and in the 
AdS/CFT correspondence [TTJI122] . 

One important consequence of the CME is the existence of a new type of collective 
propagating excitations of chiral charges: the Chiral Magnetic Wave (CMW) (221121] • In 
heavy-ion collisions, the CMW leads to a buildup of net electric quadrupolc moment in the 
plasma fireball (25H27], resulting in a charge-dependent elliptic flow of pions (2S1 EH I2H]- 
There are experimental results which appear to be in line with this picture [291T32] . 

The aim of this paper is to study these effects in a microscopic framework to gain fur¬ 
ther insight on the physics of chiral anomaly in high temperature and/or density regime. 
Although in the low-frequency regime the effects of the chiral anomaly, due to their topolog¬ 
ical nature, should not depend on the microscopic nature of the interactions, these effects in 
high-frequency regime depend on the microscopic dynamics. In this work, we study chiral 
transport phenomena in a weakly-coupled regime using kinetic theory. 

There has been recent progress on the kinetic theory with massless chiral fermions [221 
[28] . and in particular, the non-equilibrium kinetic equation describing the motion of chiral 
particles was derived. A crucial ingredient in the chiral kinetic theory is the Berry’s phase 
in momentum space which originates from the underlying chiral spinor wave function of 
quasi-particles [88] Oj. Remarkably, the effects of quantum chiral anomaly are accounted 
by the Berry’s phase in such a (semi-)classical kinetic theory. Of particular importance is 
the modification of the velocity (i.e., its relation to the momentum) of a chiral particle in 
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the presence of magnetic field due to the Berry’s phase. Such anomalous velocity will indeed 
play a central role in this paper. 

We begin by studying the collective waves in a plasma of a single Weyl fermion species in 
the regime of high temperature (and low chemical potential) in the presence of a magnetic 
held. As the collective excitations such as the CMW are hydrodynamic modes, we need to 
introduce a relaxation dynamics provided by the collision term in the Boltzmann equation, 
to be able to derive such modes in the kinetic theory. We show that the group velocity 
of the chiral magnetic wave at long wavelengths is given universally by the magnetic held 
and charge susceptibility, without depending on explicit details of the collision term. On 
the other hand, kinetic theory as a microscopic theory allows us to study the dispersion 
relation beyond hydrodynamic regime, and for this purpose, we use the relaxation time 
approximation (RTA) as a simplest example for the collision term, which enables us to find 
some analytic solutions. 

QCD and similar vector-like theories contain chiral fermions (quarks) with opposite chi¬ 
ralities - a Dirac fermion consists of a left-handed and a right-handed Weyl fermion. In non- 
abclian gauge theories, there are topological sphaleron huctuations that induce transitions 
between opposite chirality massless fermions. These non-perturbative effects, in addition to 
small quark masses, provide the dynamics of axial charge relaxation and we will refer to 
these generically as inter-chiral transitions. 

We study how the inter-chiral transitions affect the collective modes of a high temperature 
plasma in a magnetic held which contains both left- and right-handed chiral fermions. At 
sufficiently low frequency (wave number) where the inter-chiral transitions play a significant 
role, we find an interesting diffusive mode whose diffusion coefficient is governed by dynamics 
of chiral anomaly. With increasing wave number, we observe a transition of this mode into 
the chiral magnetic wave. A similar result was also found recently in a holographic model 
study [33] • 

We also explore the collective modes in the opposite regime of a cold dense Fermi liquid of 
chiral fermions [33) [35], which could be of relevance to dense quark matter in compact stars. 
In normal Fermi liquids, it is well-known that in the collisionless high-frequency regime there 
exist collective propagating modes of Fermi surface fluctuations - the Landau’s zero sound. 
The zero sound in a Fermi liquid of chiral fermions in the presence of magnetic held was first 
studied in Ref. [30] • However, only the collisionless regime neglecting collision terms in the 
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Boltzmann equation was considered. Since the CMW is a collective hydrodynamic mode 
which can be seen only with a collision term, we extend the analysis of Ref. [4Dj by including 
a collision term in the relaxation time approximation. As a result we find an interesting 
transition of the chiral magnetic wave in the hydrodynamic regime into the zero sound in 
the collisionless regime as one increases the wave number. 

The paper is organized as follows. In the next section, we give a brief review of the chiral 
kinetic theory recently developed in Ref. [331 [34] that most of our analysis is based on. In 


the subsequence section III, we study the collective modes with magnetic held in the hot 


plasma of a single right-handed Weyl fermion species. In section IV, we study the plasma of 
chiral fermions with opposite chiralities, including inter-chiral (sphaleron) transitions in the 
collision term. In section [VJ the collective waves with magnetic held in a chiral cold dense 


Fermi liquid are studied. We conclude in section VI 


II. REVIEW OF CHIRAL KINETIC THEORY 

This section gives a brief review on the chiral kinetic theory developed in Ref. [331 (3lJ, 
which also serves to set up the notations and conventions we use throughout the paper. 

In general, kinetic theory describes the motion of weakly interacting quasi-particles in 
the regime where collisions are rare enough so that the average distance traveled by quasi¬ 
particles between each collision is larger than the quantum wavelength of the particles, 
and the motions between collisions can be treated classically. This allows one to introduce 
a distribution function of quasi-particles in (x,p) phase space, and the dynamics of this 
distribution function is governed by the appropriate Boltzmann equation. We are interested 
in the kinetic theory of chiral fermions obtained from quantizing relativistic Weyl fermion 
fields. For simplicity, we restrict our discussion to the case of a single Weyl fermion held 
in this and the next section |III[ The generalization to the case of two Weyl fermions with 
opposite chiralities (or equivalently, a single massless Dirac fermion) is straightforward and 
will be considered in section IIV1 

Quantizing a single right-handed Weyl fermion held, one obtains right-handed chiral 
particle and left-handed anti-particle. With the assignment of a 17(1) charge +1 to the 
Weyl fermion held, the particles (anti-particles) carry the charge +1 (—1). We let these 
charges couple to a 17(1) gauge held (which we will call “electromagnetism”) with a weak 
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coupling strength e, which we assume to be non-dynamical 1 . Let us denote the distribution 
functions of particles and anti-particles by f±{t,x,p) where ± means particle and anti¬ 
particle respectively. The Boltzmann kinetic equation for f±(t,x,p ) reads as: 


df± ,.df± df± r ! 

x ■ —~ + P± ■ uyr = C ±l /+> /- ] » 


(II.1) 


dt dx dp 

where the equations of motion of the chiral particle and anti-particle, expressing (x,p±) 
in terms of (x, p), can be obtained from a semi-classical one-particle action derived from 
path integral formulation of chiral fermions as shown in Ref. [33]. The crucial ingredient in 
the chiral kinetic theory is the Berry’s phase and the corresponding geometric connection in 
momentum space arising from the projection of momentum dependent chiral spinor wave- 
function. In the presence of a homogeneous magnetic field B (without turning on an electric 
field), the resulting equations of motion can be written as 


VG x = v + (v ■ b)B , 
a JG p± = ±v x B , 


where 


b = 


P 

2|p| 


with 


P = 


P 

\p\ 


is the curvature of the Berry’s connection in momentum space, 

VG = 1 + B ■ b, 


(II.2a) 
(II.2b) 

(II.3) 

(H-4) 


is the modified phase space density due to the interplay between the magnetic field and the 
Berry’s curvature, and 


v = 


ck 

dp ’ 


(II.5) 


is the quasi-particle velocity with dispersion relation e = e(p). The equations of motion for 


particle and anti-particle are related by b —y —b and B —> —B, so that in Eq. (II.2a), x 
and a /G are identical for particle and anti-particle. 

In the presence of magnetic field, the Lorentz invariance dictates the existence of a spin- 
magnetic moment coupling in the dispersion relation 


e{p) = \p\ 


B ■ p 

2\p\ 2 


(II.6) 


This assumption can be justified in a theoretical limit where e —> 0 while eB kept fixed. 
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which is essential for reproducing the correct value of chiral vortical effect at weak cou¬ 
pling [SH] . 

The chiral kinetic theory is well justified only in the regime where the quasi-particle 
momentum \p\ is not much smaller than the temperature, so that the corrections to the 
dispersion relation from thermal effects can be neglected. In addition, we need to have 
\p\ \J~B for a valid classical quasi-particle picture of motion in the magnetic field, instead 
of quantized Landau level picture. This can also be seen in Eq. (II.4), where a [G as a 
modified phase space density will not make sense when y/G becomes negative in the region 
of momentum space \p\ < y/B. In the region \p\ < y/B, quantum description is needed to 
describe the physics properly. Therefore our results based on chiral kinetic theory receive 
corrections from this quantum region. Nonetheless the size and importance of the quantum 
region is suppressed for B <T 2 and the dominant contributions to the anomalous transport 
phenomena (linear in B) come from the thermal momentum region \p\ ~ T. The corrections 
from the quantum regions are of higher order in small B expansion, i.e., we expect our results 
to be valid only up to linear order in B while the quantum region contribution is of order 
B 2 or higher. We emphasize that chiral kinetic theory allows us to study a frequency and 
wave number region beyond the hydrodynamic regime, which is interesting even if we are 
restricted to the leading effects in small B expansion. 


III. CHIRAL MAGNETIC WAVE IN HOT WEYL GAS 

In this section, we consider a high temperature plasma of a single species of right-handed 
Weyl fermion (that contains right-handed particles and left-handed anti-particles) in the 
presence of a homogeneous magnetic field B to study the collective propagating mode of 
chiral charges originating from triangle anomaly: the chiral magnetic wave. We assume that 
the temperature is high enough, T /B and T>/i, such that our results are well justified 
up to linear order in B/T 2 as discussed in the previous section. This also justifies our weak 
coupling picture in a QCD like gauge theory, although we will not be specific about the 
interactions in our relaxation time approximation. In section [VJ we will discuss the opposite 
limit /r > T, of cold dense system of chiral fermions. 

For simplicity, we shall consider the case of neutral plasma p = 0 where the charge 
fluctuation modes of interest to us do not mix with the energy-momentum fluctuations. 


6 



We solve the kinetic equations Eq. (II. 1) linearized in fluctuations of f±(t,x,p ) from their 
equilibrium values, 


f°(e) = 


+ 1 ’ 

for both particle and anti-particle. Parameterizing the fluctuations with h±(t,x,p ) as 


(HI-7) 


f±(t,x,p) = / (e)- g^-h ±(t,*,p), 


(III-8) 


and linearizing the collision term (noting that the collision term vanishes for equilibrium 
distributions, C[/°,/°] = 0), 


c±[f + ,f-] = 


df(e) 

de 


l ± [h + ,h_] + 0(h 2 ), 


(III.9) 


with a functional linear operator T± acting on h±, the linearized kinetic equation takes the 
following form in Fourier space, 


iu>h± + x ■ ( ik ± B x — ) h± = X±[h + , h- ], 


(III.10) 


with frequency-momentum (cu,fc). 

In the next subsection III A[ we will derive the dispersion relation of chiral magnetic 
wave in the hydrodynamic regime of low frequency and wave number limit in a general way 
without using any explicit form of the collision operator X ± demonstrating the generality 


of the chiral magnetic wave. In the subsection III B, we will consider the case of relaxation 


time approximation for the collision term that allows us to find some analytic solutions for 
finite frequency and wave number outside of the hydrodynamic regime. 

A. Dispersion relation in hydrodynamic regime 


Since the equilibrium distribution with an arbitrary chemical potential //, 

/£(«) = 1 


e /3(€=F/d 4- 1 ’ 


(uni) 


should be a solution to the kinetic equation Eq. (II. 1), the linearized kinetic equation 


Eq. (III.10) must have a static, homogeneous (k — u — 0) zero mode solution corresponding 


to an infinitesimal shift of chemical potential. With the definition of h± in Eq. (III.8), this 


zero mode becomes h + = —h_ = const (p independent), that is, 2±[1,— 1] = 0 should be 
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satisfied. At a small but finite k, the solution h± to Eq. (III. 10) describes a hydrodynamic 


mode, whose dispersion relation will be determined in the following way. 

For future convenience, we define a linear functional, (...), which can be thought of as 
averaging over the momentum space: 


where 


and 




d 3 p 


.' dn 


J (27T) 3 ’ 

06 


T 2 

= ~ 6 + ° 


B 2 

p2 


(III.12) 


(III.13) 


(III.14) 


is the charge susceptibility at zero chemical potential. Note the normalization (1) = 1. 


Let us apply the operation of “average” defined in Eq. (III. 12) to the equation Eq. (III. 10). 


The charge conservation dictates that the collision term must satisfy 

[ VG — o, 

J p 

for any distributions f±, which implies 

((Z + -Z-)[h + ,h_]> = 0, 


(III.15) 


(III.16) 


for an arbitrary fluctuation h±. It is easy to show that the Lorentz force term in Eq. (III.10) 


vanishes upon averaging using the equation (II.2a) for x and integrating by parts, so that 
we have 

— iu>(h) + ik ■ ( xh ) = 0 , (III.17) 

after subtracting the two equations, where h = h + — /r_. In general, we need to know h up 
to normalization to find the dispersion relation from the above equation. 

However, since we know that h = const when k = 0 (the zero mode that is discussed 
above), we expect h — 1 + 0(k), and inserting this into the above gives 


oj — k ■ (x) + 0{k 2 ). 


(III.18) 


Using Eq. (II.2a) for x to compute (x), we see that the first term does not contribute as it 
is an integration of a total derivative, and the second term gives 


(III.19) 














where we have used V p • b = 2 tt5 I ' 3 \p) and /°(0) = 1/2. This finally gives the dispersion 
relation lu = v B ■ k + 0{k 2 ) with the group velocity 

duj B 


v B = 


(III.20) 


dk 4n 2 x 

In contrast with the case of sound waves whose dispersion relation is ca = c s \k\ which is 
non-analytic in k, the above chiral magnetic wave has a well-defined group velocity v B — 
(du)(k)/dk)\k = Q at zero momentum, confirming the result in Ref.[[23]. Remarkably, the group 
velocity of chiral magnetic wave v B obtained from chiral kinetic theory has no dependence 
on the details of e(p) or the collision operator X. Like sound waves, it is given solely by a 
thermodynamic property (the charge susceptibility) of the system 


B. Relaxation time approximation 


In the previous subsection, we have determined the dispersion relation of chiral magnetic 
wave from chiral kinetic theory in the hydrodynamic regime. However, kinetic theory is a 
microscopic framework and contains more information than hydrodynamics. That micro¬ 
scopic information is encoded in the details of the collision operator Z±, which eventually 
determines the dispersion relation beyond hydrodynamic regime. To see qualitative features 
of such dispersion relation beyond hydrodynamic regime with some analytic control, let us 
consider a simple example of collision operator in the relaxation time approximation (RTA) 
(see also Ref. m for another application of RTA in chiral kinetic theory): 

Z±[/i+,/i_] = —- (h± Th) (III.21) 

r v 


with a relaxation time r, and 

h= l -{h + -h.), (III.22) 

from the charge conservation constraint (see ( III.16| )). h may be interpreted as a local 
fluctuation of chemical potential. 

To simplify our analysis, let us consider only the case where k is parallel to B and 
define k via k = kB/\B\. In this case, h± is independent of the azimuthal angle in the 


transverse plane perpendicular to B, and the Lorentz force term in Eq. (III. 10) vanishes. 


Then, Eq. (III. 10) becomes 


h 


(-iuj + ix-k + T 1 )h± = ± — 
\ / r 


(III.23) 
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For real k as long as Im(u;) 7^ — r 1 , we can solve (III.23) to obtain 

--1 


h+ = ±- 


T 


-ioj + ix ■ k + T -1 


h. 


(III.24) 


The mode oj = x ■ k — ir 1 corresponds to quasiparticle excitations. We, however, are 


interested in the collective inodes which are given by Eq. (III.24). Substituting Eq. (III.24) 


into Eq. (jIII.22|), we get a self-consistent integral equation that can be used to determine 

(III.25) 


the dispersion relation, 


A = 


T 


-1 


—iu + i x ■ k + r 


-1 


= 1 


Recall the definition of (...) in (III.12). 

We will solve the above integral equation up to linear order in B that our kinetic theory 
is justified. Expanding oj(k) in B we can write 2 


— Wnorm(^) T ^anom(^) T 


(III.26) 


where oj norm (k) = O(B 0 ) and uj an om(k) = 0(B 1 ). At B = 0 Eq. (III.25) for ca norm (fc) is 
easily solved giving 3 


Wnorm(^) = — i ( t 1 — kcot(kr)) . (III.27) 

Indeed, u norm (k ) is the known dispersion relation of charge diffusion under the relaxation 
time approximation. For kr 1, cu norm (fc) can be expanded as 


Wnorm(/c) « ^ (fc r ) 2 + ^ (fc t ) 4 + • ■ • )t 4 , 


(III.28) 


and the first term gives the normal diffusion constant under the relaxation time approxima¬ 
tion, 


D = -r 

norm ^ > • 


The equation for uj anom (k) can be also solved analytically giving 

Wanom (k) = -^~k = V B k , 

4 


(III.29) 


(III.30) 


2 Since ui is parity even while B is an axial vector, the term w an0 mR') breaks parity, and is therefore labeled 
as “anomalous”. 

3 It should be noted that the solution to A = 1 exists only for kr < tt/2 because the function A(ui,k) in 


Eq. (III.251 maps the whole complex plane of to into a vertical strip |ReA| < 7r/(2 kr). For kr > 7r/2, the 


response of the gas is determined by quasiparticle excitations. 
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confirming existence of the chiral magnetic wave in the hydrodynamic regime. In the RTA 
this result for the anomalous part of the dispersion relation turns out (accidentally) to 
be exact beyond the hydrodynamic regime. More importantly, 0 J anom (k) is purely real and 
hence non-dissipative. This should remain true also outside the RTA, i.e., for a more realistic 
collision operator. 

IV. PLASMA OF FERMIONS WITH TWO OPPOSITE CHIRALITIES 

Generalizing the previous section, let us consider the more realistic case of two Weyl 
fermions of opposite chiralities (equivalently, one massless Dirac fermion), which is an ap¬ 
proximation to the QCD with light quarks. In QCD, topological sphaleron transitions can 
change the axial charge defined as a difference between the charges of the left-handed Weyl 
fermion and the right-handed one. This non-perturbative sphaleron dynamics leads relax¬ 
ation of the axial charge imbalance, whose time scale is in general different from the time 
scale of the relaxation that achieves local thermal equilibrium at fixed chirality. With the 
weak coupling QCD results in mind, we will assume that the inter-chiral relaxation time 
(which we denote as r s ) provided by sphaleron transitions is much longer than the thermal 
relaxation time r: t s 3> r. 

Based on the separation of these two scales, we can expect the following picture: On the 
time scale longer than the relaxation time r s , sphaleron transitions are effective and the axial 
charge conservation is heavily violated. What remains conserved is the total vector charge 
of the system defined as a sum of left-handed and right-handed chiral charges. We expect 
to have a simple diffusion mode of this vector charge as usual. However, at time scales 
shorter than the sphaleron relaxation time r s , but still longer than the thermal relaxation 
time r, the axial symmetry, and hence the full chiral symmetry U(1)l x U( 1)r, remains 
approximately conserved, and we expect to have the triangle anomaly in action to give us 
the features of chiral magnetic waves. There would be two propagating modes along the 
direction of magnetic held moving in opposite directions, corresponding to chiral magnetic 
waves of left-handed and right-handed Weyl fermions. In summary, as we vary the time 
scale (i.e., frequency or wave number) of the probe, we expect to see a transition from a 
simple diffusive behavior to the chiral magnetic wave. We will verify such a picture in this 
section explicitly in the relaxation time approximation. 
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The chiral kinetic equation in the case of two Weyl fermions with opposite chiralities is 


a simple generalization of our previous kinetic equation Eq. (II.1): 
f d . d . <9 \ 

\ dt + Xx ' dx + Pxi± ' dp ) /x>± = C x,±[/xT+, fx',~ } > 
where y, x' — R,L label the chirality and ± denotes the vector 1/(1) charge (particle or 


(IV.31) 


anti-particle) as before. The equations of motion Eq. (II.2) in the presence of magnetic field 
can also be generalized accordingly to 


^/G~ x x x = v x + (v x -b x )B, 


\JG X p X) ± ± v x x B , 


where we define b r = b,bi = —b, 


\[Gx — 1 + B ■ b x , 


(IV.32a) 
(IV.32b) 

(IV.33) 


and v x = de x /dp with e x = \p\ — ( B ■ b x ) \p\. The above generalization can easily be deduced 
by simple flipping of the sign of the Berry curvature for opposite chirality of fermions. 
As before, we linearize the above kinetic equation by expanding around the equilibrium 
distribution /° = l/(e^ e * + 1 ), 


fxAP^G) = f\ 


o d -f. 


de . 


hx,±(p,x,t) , 


(IV.34) 


and introducing the linearized collision operators 

df° 

V.±[ /»■,+. fx',- ] = - af X*A V,+. V,-1 + o(ft 2 ). 


(IV.35) 


We consider the relaxation time approximation for the collision operators I x ,±, assuming 
that the two relaxation dynamics corresponding to thermal equilibrium within the same 
chirality and inter-chiral sphaleron transitions are independent of each other. The former 
dynamics will introduce a linearized collision operator that is similar to the one discussed 
in the previous section, 


rthermal r 




Tiv-1 = 


^y,± -F 


T 


with 


A — 2 ^ 


A - )x > 


(IV.36) 


(IV.37) 
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where the average operation (...) is defined for each chirality as before (see (111.12)). The 
h x correspond to local fluctuations of chiral chemical potentials. On the other hand, the 
inter-chiral transition dynamics will induce a collision operator 


rinter—chiral 

L x,± 


with 


[ h x \+, h x >- ] — {h x ,± -F h) ) 

•S 


h = - ( h R + Hl ) , 


(IV.38) 


(IV.39) 


which corresponds to a local fluctuation of vector chemical potential. 

Restricting ourselves to the case of k parallel to B as before (k = kB/\B \), the linearized 
kinetic equation with the above collision operators read as 


[-iuj + ix R -k + r 1 ] h R<± = ±[t l h R — r s 1 h 5 ] , 

[-iu + ix L ■ k + r _1 ] h L> ± = ± [r _1 h L + r s _1 h 5 ] , 
where we defined h 5 as 

h^ l -{h R -h L ) , 


(IV.40a) 
(IV.40b) 


(IV.41) 


which corresponds to a fluctuation of the axial chemical potential, and introduced r = 
(r _1 + t~ 1 ) — 1 • Since we consider r s r, we have r ~ r. Solving for h Xt ± we find 

_-i 


h R .± — 


T 


—ioo + ix R ■ k + t 1 

_-i 


h 


R 


^ ) h 5 I , 


hL,± = ± 


r 


[ — iu + ix L ■ k + r _1 ] 


+ ( — I h, 


To 


(IV.42a) 
(IV.42b) 


Substituting the above into Eq. (IV.37), we obtain a coupled matrix equation 


(1 - A) -A 5 (1 - (t/t s )) 
-A 1-A(1-(t/ts)) 

where A = l(A R + Al) and A 5 = ^( A R — A R ) with 



= 0 , 


(IV.43) 


As - 

x —iu: + ix Y • fc + r _1 


(IV.44) 


w • / x 

The dispersion relation is then determined by zeros of the determinant of the matrix in 


Eq. (IV.43), which is a quadratic equation involving A and A$. From this, we will obtain 


two branches of dispersion relations, uo±{k). 
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Before presenting our numerical results, it is useful to examine various limits where 
analytic solution is possible. We first note that in the infinite r s limit, i.e. with no inter- 
chiral transitions, the equation for vanishing determinant becomes 


(1 — Ar)(1 — A l ) — 0 , 


(IV.45) 


which gives two independent branches of solutions: A R = 1 and Al = 1. Recalling that 
the dispersion relation for the case of a single right-handed Weyl fermion was obtained from 


A R = 1 in the previous section (see (III.25)), we see that this limit correctly reproduces two 
separate chiral magnetic waves from each chirality. 

Let us now consider a finite t s with u,v B k <C r _1 (recall that we already consider 
t s -1 <C t -1 ). From the equation (see (III. 19|) ), 


(x R -k) R = v B k, {x L ■ k) L = -v B k, 


(IV.46) 


we have the expansion 


A « 1 + iojT + O ((wr) 2 , (kr) 2 ) , A 5 m —iv B kr + O ((wr) 2 , (kr) 2 ) , (IV.47) 


from which the equation for the dispersion relation becomes 


w±(k) = —i 


. 1 ± y/l - (2 V B kT s ) 2 


2 T s 


(IV.48) 


There are two qualitatively different regimes depending on whether kv B r s <C 1 or kv B r s 
1. The former is the regime where the inter-chiral transitions play an important role so that 
axial charge conservation is heavily violated, while the latter is the case where one can 
neglect inter-chiral transitions and the two chiral magnetic waves with opposite chirality 


can be found. To see these explicitly, in the first case, kv B r s <C 1, expanding Eq. (IV.48) 


gives one solution u + {k) = + 0(k 2 ), which is not a hydrodynamics mode: this mode 

simply corresponds to a relaxation of axial charge due to inter-chiral transitions. The other 
solution is 

cu_(fc) = - iv 2 B T s k 2 -\ -= -iD anom k 2 -, (IV.49) 

that features a diffusion behavior with a diffusion constant 


D 


anom ^‘s i 


(IV.50) 


which corresponds to the dynamics of the conserved vector charge of the system. 
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We note that the diffusion constant for the vector charge in the presence of magnetic field 
is not determined by the thermal relaxation time t, but by the dynamics of chiral anomaly. 
This can be understood as follows. Via the Einstein relation, one may relate D anom to the 
electric conductivity 

C^anom = XU-^anom = j (IV.51) 


where xv is the vector charge susceptibility which is xv — Xr + Xl = T 2 / 3. To understand 
this, consider a homogeneous electric held E applied parallel to the magnetic held. Due to 
chiral anomaly, there is a constant generation of axial charge n§ with the rate l/(2n 2 )E • B , 
and the time evolution equation of n 5 , 


77-5 =- 77-5 + 


T, 


2n 2 


(IV.52) 


gives the stationary value 7i 5 = t s EB/( 2n 2 ) or, equivalently, the axial chemical potential 
/i 5 = t s EB /(2n 2 Xv)- With the magnetic held, this induces the chiral magnetic current 

t k B 2 


T - R - 


-E = XvV 2 B T s E = or. 


,E. 


(IV.53) 


(2tt 2 ) 2 xv 

which reproduces the above conductivity (we used v B = -B/(47 t 2 x) = B / (2n 2 xv))- Indeed, 
Eq. (IV.51) is in agreement with the result of Ref. [42] which obtained the anomalous con¬ 
ductivity cr anom in a zero temperature, finite density system using the linear response theory, 


which also has been checked in holography 

Alternatively, D anom may also be interpreted as a diffusion constant arising from the 
collective motion (drift) of chiral charges by chiral magnetic wave of the speed v B with a 
relaxation rate 1 . 

On the other hand, in the opposite regime kv B r s 1, it is easy to see that inter-chiral 


transitions do not play a role and Eq. (IV.48) produces two chiral magnetic waves moving 
in opposite directions. We see these features in our numerical solutions presented below. 

In FIG. [lj we plot the dispersion relations u±{k) determined from Eq. (IV.43) with 
v B = 0.1, t/t s = 0(red), t/t s = 0.05(blue) and t/t s = O.l(green), which feature the above 
discussed transition from the diffusive mode in kv B r s -C 1 to the chiral magnetic waves in 
kv B r s 1 . As we have anticipated earlier, for small k, both uj± are purely imaginary, and 
u: + (k ) is a non-hydrodynamic mode with a finite imaginary part even at k — 0, and c U-(k) 
is a diffusive mode with the diffusion constant given by D anom + D norm where D norm — r/3 
((III.29)) (see the right plot). Once k reaches a critical value, k ~ (2r^r s )^ 1 according 
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Re (oj r) 


- t/t s =0 


-Im(w t) 




k r 


FIG. 1: Left: Rea;_(fc) obtained from Eq. (IV.43) at r/r s = 0 (red), t/t s = 0.05 (blue) and 


t/t s = 0.1 (green). Right: -ImtJ_(fc) obtained from Eq. (IV.43) with the same values of t/t s . In 
both plots, the speed of chiral magnetic wave is fixed at vb = 0.1. 


to Eq. (IV.48), the two modes u)±(k ) meet with each other, and then split again into two 


modes now having non-zero real part (see the left figure), and each mode starts to have a 
non-zero group velocity featuring the chiral magnetic waves moving in opposite directions. 
This group velocity becomes close to the asymptotic chiral magnetic wave velocity vb as 
k increases. The (negative) imaginary part of u>±(k) are well described by D noim k 2 as k 
increases. 

We note that these features are also found in a recent work in AdS/CFT correspon¬ 
dence pSj where the authors introduced axial charge breaking (corresponding to inter-chiral 
transitions we consider) via a mass term for the axial gauge held in the holographic 5 di¬ 
mensional space-time. 


V. COLLECTIVE WAVES IN COLD CHIRAL FERMI LIQUID 

In this section, we consider a different temperature/density regime: a cold dense system 
of chiral Fermions with p T and fiT B. The second condition is needed to ensure that 
the Landau quantization can be neglected within the thin shell of width T around the Fermi 
surface where the quasi-particles can scatter. The chiral kinetic theory based on weakly 
interacting quasi-particles around the Fermi surface has been discussed recently [531 35], 
and we would like to study some of the collective modes of the system in the presence of 
the magnetic held. In particular, we show the existence of the chiral magnetic wave in this 
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system. Moreover, in the interacting cold Fermi liquid, there is another type of collective 
propagating mode - the Landau’s zero sound. We will show how the dispersion relation of 
the zero sound mode is modified by the chiral anomaly and study the transition from the 
chiral magnetic wave in the hydrodynamic regime to the zero sound mode in the collisionless 
regime as one increases the wave number k. 

At low temperature, we only need to consider quasi-particle excitations near the Fermi 
surface, and for simplicity we will restrict our discussions to the case of a single right-handed 
Weyl fermion. The chiral kinetic equation for these quasi-particles takes a similar form to 
the previous one in Eq. (II. 1) (c.f. Ref. [33]), 


dn . dn . dn , 

m +x -8i +p -d-p =c[n] ’ 


(V.54) 


where the distribution function is now denoted by n according to the prevailing convention, 
and the equations of motion for (x,p), in the absence of external electric held, are given by 


VG x = v + B (v ■ b) + £ x b , (V.55a) 

VGp = £ + vxB + (E-B)b, (V.55b) 

where \[G = 1 + B b, v = de/dp , £ = —de/dx and e = e(p) is the energy of a quasi-particle 
at momentum p. Importantly, the quasi-particle energy e(p) is also a functional of n due to 
interparticle interactions. In particular, the interaction is responsible for the effective force 
£ driven by inhomogeneity of the particle density n: £ = —(8e/5n) o (dn/dx). 


We now linearize Eq. (V.54) in fluctuations of n around the equilibrium Fermi-Dirac 
distribution 

n « = e a(.o-t) +1 ' < v - 56 > 

with the equilibrium dispersion relation e 0 = e 0 (; p)- We expand n(t,x,p), and functionals 
C[n] and e = e(p)[n], dehning function h(t,x,p) 

dn 0 . 


n = n 0 + 




(V.57) 


linearized collision functional X[h] (noting that collision integral vanishes in equilibrium 
C[n 0 ] = 0) 

C[n] = 0 + ^ X[h] + 0(h 2 )- (V.58) 
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and linearized Landau interaction functional J-[h\ 


e(p)[n\ = e 0 {p) + F[h\ + 0(h 2 


(V.59) 


The conventional choice of the definition of the function h in Eq. (V.57) is such that h = const 
describes a change of the distribution due to a shift of the chemical potential /i by h. We 
also note that in the zero-temperature limit (j3 —» oo) the coefficient dn 0 /dn —> <5(eo — /i), 
i.e., the fluctuations of n are constrained to a thin shell around the Fermi surface. 

As discussed in Ref. [35] from the point of view of quantum field theory, Eq. (V.54) and 


Eq. (V.55) may be derived from high density effective theory in which the chemical potential 
fi is much larger than the excitation energy of quasi-particles of order T and a J~B. Since the 
derivation is only up to a linear order in B, we will keep our results only up to this order 
for consistency. 


Expanding the second term on the l.h.s. of the kinetic equation (V.54) using equation of 

dh dn$ 


motion (V.55a) we find 


. dn 

X • 7 — = Xo ■ 

dx dx 


+ 0(h 


where we denote by £ 0 full quasiparticle velocity in equilibrium 

1 


x 0 = 


Vg 


(«o + B(v 0 • 6 )), 


(V.60) 


(V.61) 


with Vq = deo/dp being the normal particle velocity in equilibrium. We linearize the third 


term on the l.h.s. of Eq. (V.54) using equation of motion (V.55b) and definitions (V.57), 


(V.59) and (V.61) 


P 


dn 

dp 


-£ ■ x 0 + (x {) x B ) • ^(1 + F)[h} 


dno 

dp 


+ 0(li 2 


(V.62) 


We used the fact that £ = 0{h ) since, by definition, £ = — dF[h]/dx + 0{h 2 ) and applied 
the relation dno/dp = —Vo(dn 0 /dn)- 

Substituting Eqs. ( V.60[ ) and (V.62) into Eq. (V.54) one arrives at the linearized kinetic 
equation, which after Fourier transform from x to k takes the form 


iojh + x 0 • ( ik + B x — ) (1 + F)[h] = Z[h\. 


(V. 63) 


A. CMW at small wave number 


Since the Fermi-Dirac distribution no in Eq. (V.56) is a solution to the kinetic equation 


(V.54) for any /r, the linearized kinetic equation (V.63) must have a zero mode at k = lu = 0 
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In the noninteracting Fermi gas, that mode would be simply h = const (p- independent), 


since h, as defined in Eq. (V.57) coincides with the shift of the chemical potential. Interac¬ 


tions, however, cause the change of the energy €o(p) to be different as well, and thus the two 
distributions at chemical potentials /i and p! have different functions eo (p) and e' 0 (p) in the 


Fermi-Dirac distribution Eq. (V.56). The difference between two such infinitesimally close 


distributions defines Sh - the zero mode, according to Eq. (V.57), 

dn 0 


<9/i 


Sh = n 0 - n 0 = ^(<5/1 - <5e 0 ) • 


(V.64) 


Since Se o = e' 0 — eo = J-[5h], we find that even if Sh could depend on p, the functional 


Sh + ^[Sh] must be p-independent, or 


(1 + SF)[ho\ — 1, 


(V.65) 


where we defined (normalized) the zero mode as the ratio of two infinitesimal quantities: 
ho = Sh/Sp. In the absence of interaction, T = 0, the zero mode is ho = 1. 


From the fact that C[n' 0 ] = C[rio\ = 0 and the definition of X in Eq. (V.58) it follows that 


X[ho\ — 0. 


(V.66) 


The left hand side of the kinetic equation (V.63) also vanishes for h = h 0 due to Eq. (V.65) 
at u> = k = 0. 


At small but non-zero k, the solution h to Eq. (V.63) describes a hydrodynamic mode, 


whose dispersion relation we can determine in the following way. 

For further convenience we define a linear functional, (...), which can be thought of as 
averaging over the momentum space 

<9 no 


XoJp <9/i 


where 


Xo 


= [Vg 

J n 


dp 


(V.67) 


is the susceptibility of a free Fermi gas (neglecting integraction). 


Let us apply operation of “averaging” defined in Eq.(V.67) to the equation (V.63). Con¬ 


servation of charge means f \/~GC\n ] = 0 for any n and therefore requires that (X[/i]) = 0 for 


all h. The last term on the r.li.s. of Eq. (|V.63|) also vanishes upon such averaging, because 

'd n 0 




dp 


xB(l+J)[ft] =0 (V.68) 
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is an integral of a total derivative. Thus we obtain, for any h, 


— iu(h) + ik • (x 0 (l + J r )[/i]) = 0. 


(V.69) 


In order to determine the dispersion relation of the mode, in general, we would need to know 
h(k,p), corresponding to this mode. However, for k = 0 we already know that the mode is 


given by ho- Using Eq. (V.65), we can then write the frequency of the hydrodynamic mode 
to leading order in k as 

/ 

(V.70) 


C0 = k .^l + O(k 2 ) 


(ho) 


Using the definition of the full velocity Xq in Eq. (V.61) we see that the first term does not 


contribute to the average ( x 0 ), since it gives an integral of a total derivative, and as a result 


... If dn 0 , _. B f 

(X 0 ) = — (To + B(v o -b)) = — / 

Xo Jp Xo Jp 


drin 


• b = 


B 


(V.71) 


dp J “ 4tt 2 xo ’ 

where we integrated by parts and used (d/dp) • b = 2tt5 3 (p) and no(p = 0) = 1. The 
coefficient 1/47r 2 here is the same as in the expression for the chiral magnetic effect J = 
1/(47 T 2 )flB. 

The velocity of the hydrodynamic mode (chiral magnetic wave) is given by 

doj B 


V CMW — 


(V.72) 


dk 47r 2 xo(h 0 ) 

The product x 0 (h 0 ) is the full susceptibility x °f the Fermi liquid with the interactions. 
Indeed, by definition, x is the rate of change of the volume density riy as a function of the 
chemical potential, 


X = 


5n 


v 


Sji 


(V.73) 


where ny = f y/Gn. Thus, using Eqs. (V.64), (V.67) and the definition of h 0 = 5h/Sp, 


X 


= Vg 


drin 5h 


= Xo(K) ■ 


(V.74) 


Jp dp 5p 

This shows that the chiral magnetic wave velocity is still given by the same expression in 
terms of susceptibility in the cold dense Weyl liquid. As in the case of finite temperature, 
it also has no dependence on the details of dispersion relation e(p). Remarkably, the same 
relation between the speed of chiral magnetic wave and the susceptibility still holds for Fermi 
liquid even if the susceptibility is modified due to interaction. 
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B. Transition of CMW into the zero sound 


In the previous subsection, we have established the chiral magnetic wave in a cold Fermi 
liquid system in the hydrodynamic regime io,k <C r _1 , where r is the relaxation time due 
to particle collisions. In the opposite limit of io,k r _1 , the Fermi liquid has another 
collective excitation - the zero sound. At low temperature, a typical relaxation time scales 
as r ~ h-/T 2 due to Pauli blocking effect [3Ji 05], and if the temperature is strictly zero, 
the relaxation time diverges and one is always in collisionless (zero sound) regime. We will 
consider a small, but finite temperature so that one can discuss the transition from the chiral 
magnetic wave at w C r' 1 to the zero sound at to r _1 . 

For this purpose, let us study the dispersion relation in the kinetic theory within the 
relaxation time approximation, 

l[h] = _ h ~ ' (V.75) 

T 

The general form of the linearized Landau interaction functional F[h\ to linear order in 
B was discussed in Ref. [3]JJ. As shown in Ref. [4E], f° r massless quark matter weakly 
interacting via exchange of color gluons in the absence of external magnetic field, one can 
approximate F[h\ = Fq(Ii ) with a constant F 0 to lowest order in a s . For simplicity, we will 
assume this form of the Landau functional in the following analysis. 

The linearized kinetic equation with these assumptions reads 

(— iio + ix o ■ k + r -1 ) h = ( t -1 — iF 0 Xq ■ k ) ( h) , (V.76) 


where the Lorentz force term disappeared as we assumed, for simplicity, k parallel to B as 
before. The equation for the dispersion relation is 

■ -1 - iF 0 x 0 ■ k 


T 


-iu + ix o • k + t~ 1 


= 1 . 


(V.77) 


It is easy to see that when r _1 —» oo (i.e., lot 1) we recover equation (V.69) obtained 
earlier in hydrodynamic limit. In the opposite, collisionless limit, lot 1, we can drop r^ 1 
and the resulting equation is the well-known equation for the zero sound. 

In solving this equation, we consider the temperature low enough that dno/dfi can be 
replaced by 5(eo — n) up to higher order corrections in T / /i. Let us first consider the case 
without magnetic field in the collisionless regime lo t _1 ; in this case it is easy to solve 
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(a) 


(b) 


FIG. 2: The numerical solution of Eq. (V.77) for the real and imaginary parts of the collective 
mode frequency u(k) as functions of the wave number (in units of t -1 ) in a cold Fermi liquid with 
Landau parameter Fq = 1 and vb = B/4tt 2 x = 0.1. The dashed curves are the same but in zero 
magnetic field. 


(V.77) to obtain 


where s 0 is a solution of 


w 0 (fc) = s 0 v F k, 


v F = 


de 0 

dp 


> 


60 =/J- 



(V.78) 


(V.79) 


This is the well-known result for the zero sound in normal Fermi liquid |44j . Since we have 
shown in small B limit the existence of chiral magnetic wave in the k <C r _1 regime, we 
naturally expect a transition between the chiral magnetic wave and the zero sound as k 
increases from k <C r _1 to k r _1 , which we confirm in our numerical results. 

In FIG. [2j we present our numerical results for the dispersion relation u>(k) with F 0 = 1 
and v F = -j=?- = 0.1. Results with a different choice of (. F 0 ,B ) are qualitatively similar. 
For comparison, we also plot the dispersion relation without the magnetic field (the dashed 
curve). As one can anticipate, oj(k) without magnetic field is purely imaginary for sufficiently 
small k due to diffusion, while for a large k one finds the propagating zero sound with 
frequency-independent attenuation. For fl ^ 0, there is a linear dispersion of the CMW for 
small k t -1 , and as k increases we observe the transition to the zero sound, modified by 
the magnetic field. 
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VI. SUMMARY AND CONCLUSION 


In this paper, we studied collective modes of a weakly coupled plasma of chiral fermion 
quasiparticles in the presence of an external magnetic field, using the framework of the 
recently developed chiral kinetic theory. We show the existence of the chiral magnetic wave 
in both the high temperature, low density regime, as well as the cold, dense Fermi liquid 
case, in a general way without relying on an explicit form of the collision term. The velocity 
of the chiral magnetic wave is shown to be given universally by vb = B/(47r 2 x) where x is 
the charge susceptibility, even when the interactions in the Landau Fermi liquid modify the 
quasi-particle energies. 

In a gas of massless Dirac fermions, we show how inter-chiral transitions cause chiral mag¬ 
netic wave to become a vector-like charge diffusion mode at sufficiently low wave numbers. 
In the interacting cold Landau Fermi liquid, we observe an interesting transition from the 
chiral magnetic wave to the zero sound as the wave number increases from the hydrodynamic 
to the collisionless regime. 

One of the important approximations in this work is the decoupling of the charge fluctua¬ 
tions we consider from the energy-momentum fluctuations. In high temperature, low density 
regime (discussed in the sections [TTT] and IV), it can be justified as far as T > /i. In high 


density, low temperature regime (discussed in the section 0, this assumption may again be 
justified for large N theories as the coupling between the charge and the energy-momentum 
fluctuations is 1/N suppressed. It would be interesting to extend our analysis including the 
energy-momentum fluctuations to see, for example, the mixing between the first sound and 
the chiral magnetic wave. Another interesting direction is to let the electromagnetism be 
dynamical, and study dynamical collective excitations such as plasmons |23U17] as well as 
their interplay with the chiral magnetic wave. We leave these and other extensions of this 
work to future studies. 
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